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We calculate lepton flavor violating (LFV) Z decays Z — > e^ej = e, fj,,T ; i ^= j) in the 
Zee model keeping in view the radiative leptonic decays et — > ej7 ( i — n, t ; j = e, n ; i 7^ j ), n 
decay and anomalous muon magnetic moment QuAMM). We investigate three different cases of Zee 
fij coupling (A) f% = f% = f 2 Te , (B) f% » f? e » fl T , and (C) f% > / 2 M > f 2 Te subject to the 
neutrino phenomenology. Interestingly, we find that, although the case (C) satisfies the large excess 
value of fiAMM, however, it is unable to explain the solar neutrino experimental result, whereas the 
case (B) satisfies the bi-maximal neutrino mixing scenario, but confronts with the result of /1AMM 
experiment. We also find that among all the three cases, only the case (C) gives rise to largest 
contribution to the ratio B(Z —* e t^) / B(Z — > ft /x^) ~ 10 -8 which is still two order less than 
the accessible value to be probed by the future linear colliders, whereas for the other two cases, this 
ratio is too low to be observed even in the near future for all possible LFV Z decay modes. 

PACS number(s): I3.38.Dg, 13.35.-r, I4.60.-z, 14.60.Pq. 



I. INTRODUCTION 

The high statistics results of the SuperKamiokande 
(SK) atmospheric neutrino experiment [jjj and the so- 
lar neutrino experiment (2| have strengthen the conjec- 
ture of neutrino flavor oscillation from one species to an- 
other. The phenomena of neutrino oscillation leads to 
non-zero neutrino mass and the scale of which is ~ eV 
predominantly set by the atmospheric and solar neutrino 
experimental results. Such a tiny neutrino mass could 
be generated by several ways, namely, see-saw mecha- 
nism H , non-renormalizable operators [Q or through the 
radiative ways at the one or two loop level. One of the 
most well known model of radiative neutrino mass gener- 
ation is proposed by Zee Q in which small neutrino mass 
is generated at the one loop level due to charged scalar 
exchange through explicit lepton number violation. The 
model has been investigated by many authors [6 - 11]. 
The model contains one extra charged singlet scalar field 
with non-zero lepton number and another doublet Higgs 
field in addition with the standard model (SM) contents. 
The scalar field content of the Zee model is not only re- 
sponsible to generate tiny neutrino masses but also gives 
rise to non-standard interactions due to the presence of 
charged scalar fields, e.g., one of them is the anomalous 
muon magnetic moment (/xAMM). The excess value of 
(/zAMM), Sa^ = (43± 16) x 1CT 10 recently reported by 
E821 Collaboration |l2] leads to the new source of in- 
teractions beyond the SM level. Furthermore, the Zee 
model leads to possible Lepton Flavor Violating (LFV) 
Z decays such as, Z — > efej (i, j = e, jtz, r and 



hereafter, we will assume i 7^ j unless otherwise stated ), 
which have taken interest in view of future collider plans. 
The present sensitivity of the measurement of branching 
ratio of Z — > efef at LEP is ~ 10 -5 whereas future lin- 
ear colliders (NLC, JLC, Tesla GigaZ) will bring it down 
to ~ 1CP 8 and thus the testability of such model will be 
increased due to higher sensitivity of measurement which 
could be able to reveal new physics beyond the SM. 



In the present work, we calculate /1AMM and LFV 
Z decays in the Zee model keeping in view the other 
constraints arising due to \x decay and other LFV ra- 
diative lepton decays. We estimate pAMM and LFV 
Z decays by utilizing the constraints on the parameter 
space obtained from the /1 — > eVeV^ decay and also from 
radiative — > ejj ( i — /1, r ; j — e, fi ; i 7^ j) de- 
cays. Recent works in this path have been done H] in 
which the Zee model have been investigated in view of 
recent /iAMM experimental result and LFV decays. In 
the present work, we restrict ourselves within the config- 
uration of the minimal Zee model and we, particularly, 
investigate different hierarchical cases of Zee /y coupling 
subject to the present neutrino phenomenology. The plan 
of the paper is as follows: In Section II, we will first briefly 
review the Zee model, its basic interaction Lagrangian, 
the charged scalar mixing and the neutrino mass matrix. 
The constraints on the Zee /y coupling due to /a — > eVeV^ 
, &i — > ej7 decays and /j,AMM are discussed in Section 
III. The LFV Z decays are calculated in Section IV and 
Section V contains summary of the present work. 
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II. BRIEF REVIEW OF THE ZEE MODEL 



A. The interaction Lagrangian 



where we have dropped the quark interaction terms, and 
lcpton number conserving Higgs potential terms. The 
lepton doublets are denoted as Iil {i = 1, 2,3) with the 



definition l c iL = (^l) c = C£ iL and <fi a (a = 1, 2) are the 
The interaction Lagrangian of the Zee model is given Higgs doublets and h~ is a charged singlet scalar field. 



by 

C = ^2 fatiLiT2t c jL h~ +c 12 (f)iiT2<j)2h~ +^y l 1 l L(t>2e i R+h.c. 

(2.1) 



The fij terms can explicitly be written as follows: 



E 



fijtiLiT2t C jL h + h.C. = ^ fiji^iL^jL ~ e iL^ L ) h + h -C 



= 1\f eli {v eLl i c L - e L v^ L ) + f^iy^rl - H l v c tL ) + f T e(^rLe c L - T L v c eL )}h + h.c. . 



(2.2) 



Since eiLV c , L = VjLC c iL , we can also re-express (2.2) as 



2[/e/x(ML^eL - e L^ L ) + Ut(tl^l ~ Vl v tl) + Jre(eL^L - T L V° eL )\h + h -C 



(2.3) 



The C12 term can explicitly be expressed as follows: 

ci2<^V 2 <M~ + h.c. = c^t '4>l - 4>\4>t)h~ + h.c. (2.4) 
We define 



(2.5) 



X° = x?cos/3- X Osin/3, x° = X? sin/3+x" cos/3 , (2.6) 

H° = flj cos/3 - Hi sin j3 , H° = H° sin/3 + ff 2 ° cos/3 , 

(2.7) 

4> + = <f>i cos [3 — <f>2 sin /3 , 4> + — <fit sin /3 + (f>2 cos /3 , 

(2.8) 



tan ,5 = vi/v2 , 



(2.9) 



where v a are the vacuum expectation values (VEV's) of 
4>° a ,(<fP a )=v a /& 

By using the expressions (2.6)-(2.9), the c\2 term is 
expressed as follows: 



c 12 (<t>i<t%-<$<fi)h- + h.c. 



ci2 [<f>+ {v + H°- ix°) -<P+{v + H°- i X °)]h- + h.c, 

(2.10) 



V2 
where 



v = v\ cos (3 — v 2 sin /3 = 0, 



(2.11) 



v = vi sin /3 + v 2 cos [3 = \/vf + v\ 



2M W r- 

= V2 x 174GeV. 
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(2.12) 



The components </> + and x° are absorbed into the gauge 
bosons W + and Z°, so that they are not physical par- 
ticles. Therefore, the physical part in the ci 2 term is 
only 

-^cuvep+h- + -^c 12 (t) + H h- + h.c. (2.13) 
V2 V2 

The first term of Eq.(2.13) is a quadratic interaction term 
which induces the <f> + -h + mixing as we show in the later 
part of this section. The second term is not relevant for 
our analysis. 

The third term in the Lagrangian (2.1) is the conven- 
tional Yukawa coupling term : 

^2vi{v iL e iL ) ( ) e iR = ^2vi( v iL4>t + c iL (j) 2 \)c iR 



2 



= -y=v cos yj (e iL e iR ) - sin /3 ^ yijviLeni)^ 

* i i 

-^X>(ei£e«)[sin/3(iT - *X°) - cos/3iJ°] 

+ (unphysical terms), (2-14) 
so that the charged lepton masses m ei are given by 

m ei = -^=yiiicosf3. (2.15) 

The term J2i Hi ^-^ipiL^iR)^ together with the terms 
given in Eqs.(2.3) and (2.13) will contribute to the neu- 
trino mass generation. 



C. The neutrino mass matrix 



The neutrino mass is generated in the Zee model due to 
the charged scalar exchange at the one loop level through 
explicit lepton number violation. The Zee neutrino mass 
matrix M v is given by the form as 

/0 p a\ 
M v =\ p A (2.23) 
\ a A / 

where p = {M v ) e/X , a = (M v ) er , A = (M„) MT and 

[M v )ij = m fcj(m e] 2 - m e 2 )/m 2 T (2.24) 



B. The charged scalar boson mixing 

The mass matrix for (h + , cf> + ) is given by 



M 2 = 



Ml c 12 v/V2 \ 
c 12 v/V2 M 2 J 



(2.16) 



where M 2 and M| are the coefficients of the h + h~ , (f) + (f>~~ 
terms of the scalar potential. We define an orthogonal 
transformation between the h + and d> + scalar fields as 



cos v sm t 
— sin cos ( 



Hi 



(2.17) 



M 2 c 12 v/V2 



by which we obtain 
= ( H+ H} ) 



(M 2 


Mia \ 




\M 12 


M 2 ) 


m 



(2.18) 



, M 2 +M 2 Ml -Ml v 
2 2 y2 



(2.19) 

_ M 2 + M 2 M 2 -M 2 y 
M 2 = h_ 1 cos 20 + ci2 sin 20, 

(2.20) 

M 12 = — £sin20 + ci2-^cos20. (2.21) 

2 V2 

The diagonalization condition gives 

2ci 2 i>/v / 2 -2\f2c 12 m w jg 



tan 20 



f * M » v/(M? - Af 2 2 ) 2 - (2V2c 12 m w /g) 2 

(2.22) 



sin 20 tan /3 m 2 Ml 
32tt 2 w/x/2 M 2 2 



(2.25) 



and sin 20 and tan/3 are determined from Eqs.(2.22) and 
(2.9), respectively. 

It has been shown in Ref.[6] that the model can ac- 
commodate the atmospheric neutrino experimental re- 
sult as well as candidature of neutrino as a hot dark 
matter component through the choice of coupling as 
fep, ~ //it > fer and the neutrino mixing matrix obtained 
as 



1 



U 









1 


1 


V2 


V2 


1 


1 


•J2 


V2 



(2.26) 



and such mixing matrix cannot accommodate solar neu- 
trino experimental results. It has been advocated to add 
a sterile neutrino in the model to explain the solar neu- 
trino experimental results ||. 

However, another interesting option to explain both 
the solar and atmospheric neutrino experimental results, 
namely, bi-maximal neutrino mixing can arise in the Zee 
model [7] due to the choice of coupling as / eM ^> 
fer 3> f i_lt which leads to the mixing matrix as 



/7f 



U 




(2.27) 



and the mass matrix resemblance to the form presented 
in Ref.[15]. An explicit realization of the above men- 
tioned hierarchy of Yukawa coupling has been demon- 
strated in Ref. [9] due to the inclusion of a badly broken 
SU(3)h horizontal symmetry through a simple ansatz 
on the symmetry breaking effects. The SU(3)h symme- 
try breaking has been considered to be proportional to 



3 



the transition matrix elements of the mass-matrix and 
thereby obtained the three Yukawa couplings 



fen = [mr/irrift + m e )]f, 
f eT = -[m^/(m T + m e )]/, 



(2.28) 



which is in excellent agreement with the experimental re- 
sults. However, in the present work, we will show that 
bi-maximal mixing hierarchy confronts with the hierar- 
chy required to explain the excess value of /iAMM in the 
Zee model. 



fur = [m e /(m T + 77i M )]/, 

which necessarily leads to the hierarchy required to ex- 
plain the solar and atmospheric neutrino experimental 
results through bi-maximal mixing. The solar and the 
atmospheric neutrino experimental results are connected 
through the relation 



Am 2 sol /Am 2 atm ~ V2m e /"V = 6-7 x 10" 



(2.29) 



III. CONSTRAINT FROM p -» eV e v„ AND e< -> e j7 
DECAYS 



A. n—>ev e Vfj, decay 



From the interactions (2.1), we obtain the effective four 
Fermi interaction 



Ceff = =5 l-2U(e L ^ L ) + 2f Te (e L ^ L )}[2f eii (V c eL f, L ) - 2f^ Ll x L )\ 
M 



=^[-/e M (ei^L)Ki/il) - frefp,T{eLVrL)(y C TL liL) 
M 



+ frefe^LVrh)(y C eL [lL) + /e M / M r (SL^L ) (Kl^L )] + 



(3.1) 



where 



1 cos 2 6 sin 2 1 



Ml 



(3.2) 



By using the Fierz transformation and the formula 

feOfe = ^(C^OCf^L, (3.3) 
we can rewrite the first term of Eq.(3.1) as 

£e// = ^=2/e M (eL7M Jy ei)( I7 Mi7' 1 ML)- (3.4) 

2 M 

On the other hand, the conventional interaction for the 
[i— *ev e Vp, decay is given by 

C = ^(e7 M (l " TsHX^a - 7s)M) + h.c. , (3.5) 



G F _ g 2 _ 1 

V2 8M 2 , 4(v/v / 2) 2 ' 



(3.6) 



so that the relative ratio of the contribution from (3.4) 
to that from (3.6) is given by 



4/ e y8M 
G F /V2 



M 



(3.7) 



The four-Fermi coupling in the Zee model comes out 
as (Gf/v / 2)(1 + C) where the parameter £ is deter- 
mined from the deviation between the observed "Gf" 
value from /i — * eF e i^ and that from hadronic weak de- 
cays. Smirnov and Tanimoto || have put the constraint 
C < 10~ 3 (so that the model does not destroy the agree- 
ment in the electroweak precision tests), and they have 
obtained 



where 



< 1.37 X 1(T 4 , 



1 M 



(3.8) 



(3.9) 



Furthermore, the second, third and the fourth terms of 
the effective Lagrangian give rise to the other possible 
/x decay modes as /i — > v T ev T , fi — > V e ev T , /i — > u T eVfj, 
respectively, which also give the same final state signal 
as ji — ► e + "missing energy" . In a similar way, the 
constraints obtained from those processes are estimated 
with C 2 < 1CT 3 as 



s/we/er; s/ur/er) s/we/ur < 4 X 10 



(3.10) 
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B. 



ej7 decay 



The non-standard radiative decay |L6[ ej — > ej7 arises 
in the Zee model due to 7-emission from the Zee boson, 



7-emission from the initial state lepton and 7-emission 
from the final state lepton. In general, the gauge invari- 
ant total amplitude of the above process is obtained as 



A(ei -> e,-7) 



e fik fkj _ / \ r, 

96ir 2 M 



rrij) + (mi - m j )^ 5 \ui(p 1 )e^(q)q 1 ' 



(3.11) 



with q = pi — P2 and the decay width as 



r(e< 



il) = 



48tt 192tt 3 V M 



fik fkj 



(3.12) 



where we have considered rnf 3> m?. To estimate the 
values of we consider the ratio of the branching ratios 
of the decays ej — > e 3 l and — > ejVjVi as 



B( ei -► ej -7) Zo 



/ fik fkj 



B{ei -> e 3 V 3 v.i) SM 6tt V [M/(w/V2)] 5 



(3.13) 



The present experimental status of the above ratio is 
given by @ 



/7 < 1.2 X 10 -11 , 



< 



2.7 x 10" 



B(t — > 67) 
B(r -> eF e j/ T ) " (17.83 ±0.06) x 10" 



< 



1.1 x 10~ 6 



5(t -> /zF M z/ T ) " (17.37 ± 0.07) x 10" 



(3.14) 
(3.15) 
(3.16) 



Therefore, we obtain the constraints on the coupling con- 
stants as 



£|/ M r/re| < 4.67 X lO" 5 , 

Z\f*M < 5.24 x 10- 2 , 



Z\Uefer\ < 3.39 X 1(T 



(3.17) 



We summarized the results obtained in Eqs.(3.8), (3.10) 
and (3.17) in Table I. 



C. Anomalous muon magnetic moment (fiAMM) 



The excess value of /iAMM recently reported by the E821 
collaboration based on the theoretical calculation pre- 
sented in Ref. [18] indicates the signal of new physics be- 
yond the SM. (However, it has been argued in Ref. [19] 
that the discrepancy between the theoretical and exper- 
imental values could be removed if other estimation of 



hadronic contribution to the photon propagator is con- 
sidered ). The excess value of 



5a, 



(43 ± 16) x 10 



-10 



(3.18) 



could arise in the Zee model due to charged scalar ex- 
change at the one loop level. The diagrams are essentially 
same as the ei — > ej^ process by regarding i — j = p. 
The extra contribution to <5a M arises in the Zee model 
due to charged scalar exchange at the one loop level and 
it is estimated as 



6a L 



^-^z— e.r f i 



2 2 



247r 2 Af| 



(3.19) 



Therefore, if we regard that the value of (3.19) comes 
from the extra contribution due to charged scalar ex- 
change we obtain the constraints on the couplings as 



ZUle +flr)= 7.67X10- 



(3.20) 
The 



where we have considered the central value of <5a M . 
previously obtained value of £/ 2 e from /i — > V^ev^ decay 
given in Eq.(3.8) is too low to explain the large value of 
Sa^, whereas the higher value of £/ 2 T ~ 7.67 x 10 _1 is not 
in conflict with the previous results given in Eqs.(3.10) 
and (3.17). It is to be noted that the hierarchy of cou- 
pling required to explain such an excess value of /iAMM 
{fur 3> fen) confronts with the hierarchy required to 
explain bi-maximal neutrino mixing as discussed in Sec- 
tion II. C. Thus, the upshot of our analysis is that if we 
consider seriously the experimental result of /iAMM ex- 
periment we have to give up bi-maximal neutrino mixing 
pattern in the Zee model. 



IV. LFV Z DECAY'S 



etej 



In this section, we calculate LFV Z decays, Z 
(i,j = e,fi,T and i ^ j) which have taken much interest 
in view of future colliders [|l3| . The explicit LFV coupling 
of the Zee model gives rise to such processes which is an 
immediate consequence of the non-zero neutrino mass. 
The diagrams of LFV Z decay process have been given 
in Fig. 1. The total contribution of all diagrams can be 
written as 



5 



A(Z(p) eiMejfa)) = ^^KuMh^l l5 )v t ( Pl )e»(p) 

Z COS Uw 



(4.1) 



with p = pi + P2 and 
K = {-2 sin 2 ^[cos 2 + sin 2 6-Ffa)]} 



and 



6 = 



M 2 



6 



+ cos 2 e-G(^) + sin 2 e-G(&) - sin 2 cos 2 6» • i2"(£i,£ 2 )} • SinCe the mixin S an S le 6 can also be written as 

(4.2) 



F(£i) = [ dx [ dy x ln[l - &zy(l - y)] , (4.3) 
Jo Jo 

Gte) = i + i ^ ^lnfl-ar-^Vl-y)] 



sin 



Mj - M 2 _ ft/fr _ l 
M 2 - M 2 &/£i - 1 



where 



[ dx [ dy- — 

Jo Jo 1 - a; - 



^x 3 y(l - y) 



£,ix 2 y{\ - y) ' 



(4.4) 



and is related to 



£ ' £0 — £1 ' £2 j 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



#(£1, £2) - ^(£1, £2) + f(£2,£i) f(£i) F(&) , (4.5) thc factor K ig g . ven ag a function of the p aramet e r s ^, 

£2 and £. The approximate expressions of F(£i), G(£i) 



^(£i,£ 2 ) = / <*r / ^ sin 
Jo Jo 



1 - (1 - -r)y -£\xy{l -y) 
£2 



and i?(£i, £2) are given in Appendix A. For Af 2 > M » 
M| > M|, a dominant term in Eq.(4.2) is thc #(£i,£ 2 ) 
(4.6) term and we get 



K ~ 7 sin 2 9 cos 2 0JT(f 1 , £2 

2(4tt) 2 



_1 frW*.^ ^_£ 

2(47t) 2 (6-£i)H£0 A £0 



0) ./) ./) 



dy xln 



1 + 



(£ 2 -ar 



2/(1 - y) 



£i£ 2 [l-£ixy(l-y)][l-6^(l-J/)] 



(4.11) 



The conventional lepton number conserving Z decay is described by 
I 



A(Z(p) -» ei(pi)ei(p2)) 



Ui(p2h»(gv - gA75)^(pi)e M (rt , 



2 cos 0jy 



(4.12) 



G 



9v = --(1 - 4sin 2 6»n/) , 3a 



r(z 



T _ G F M% , 2 



(4.13) 



(4.14) 



Therefore, the ratio of the branching ratios is given by 
B(Z ->■ - ± - =F 



Rij — 



B(Z^efef) 



2 /^x2 



(4.15) 



The experimental values of B(Z — > e T e f) are as follows 
: B(Z -> e ± /i^) < 1.7 x 10~ 6 , S(Z -> e ± r = F) < 
9.8 x 10~ 6 , S(Z -> m ± t =f ) < 1-2 x 10~ 5 . Utilizing 
these experimental limits, we obtain Rf^ p < 5.05 x 10~ D , 
R%% v < 2.91 x 10 -4 , R^ T P < 3.56 x 10" 4 by using the 
branching ratio B(Z -> f^fi*) = 3.369 x 10~ 2 0. 

Since our interest is in the maximum value of Rij un- 
der the constraints obtained in Eqs.(3.10) and (3.17), we 
have parameterized Rij in terms {£,fikfjk) 2 {K/^) 2 not by 
(fikfjk) 2 K 2 . As seen in Eq.(4.11), in the limit of £ — > 0, 
the factor K becomes vanishing while the value of \K/£\ 
increases as the value of £i decreases because K/£ is ex- 
pressed as 



K _ 

T ~ 2(4tt) 2 6 
where 
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Mi 
Ml 



(4.16) 



(4.17) 



Therefore, for a small £i (but sizable £2)1 the factor |if/£| 
is approximately proportional to m(£2/£i) an d (1 — «)■ 
Since £2 = £0/^ , the value of has a maximum 

\K/£\max at k = k — 1/2 for a large £2/^1- We illus- 
trate the behavior of |i<7£| versus k for three different 
values of M , M = 200 GeV, 500 GeV and 1 TeV (cor- 
responding £0 =0.21, 0.033 and 0.0083 , respectively) 
in Fig. 2. Furthermore, for £ > 10~ 2 , the value of kq 
deviates from ko ~ 0.5. This is due to the fact that the 
approximate expression (4.16) is valid only for £ -C 10~ 2 . 
For £ < 10 -2 , we can find that the expression (4.16) is in 
agreement with the direct numerical estimate within the 
deviation of 5 %. 

The value of \K/£\ is also highly dependent on the 
value of £. Since we want to obtain a value of \K/£\ as 
large as possible, we consider the value of £ as small as 
possible. We will estimate the lower bound of £ by us- 
ing the perturbative unitarity bound on /y coupling as 
/j 2 /47r < 1 along with the experimental constraints on 
V^fij given in Table II. It is to be noted that the exper- 
imental constraints on vf/ij depend on the assumptions 
for the hierarchical structures of fy. In Table II, we have 
considered the following three typical cases: (A) / 2 ^ = 
/ 2 T = / 2 e , (B) / 2 M » / 2 e » / 2 T , and (C) / 2 T » f% » 



/ 2 e . The case (B) is motivated by the simultaneous ex- 
planation of the solar and atmospheric neutrino data as 
we discussed in Section II. C. The case (C) is motivated 
by the explanation of the excess value of juAMM. 

In order to obtain a value of \K/£\ max at n — kq as 
large as possible, we want to take a value of M^ as large 
as possible. However, it is unlikely that the value of 
is far from the electroweak scale. For numerical estimate 

500 GeV. 

• a — ffur = fre- thi s situation, the most strin- 
gent bound comes from /1 — ► ej decay as 



of \K/£\ max , we will take M, 
(A) fl ~ 



yfl\U\ = yft\Ur\ = V^|/re| < 0.68 X 10" 



(4.18) 



We have listed those bounds on the term for all 

the three cases in Table II. In order to obtain a value 
of Rij as large as possible, we take the maximal value 
as v^l/y'l ~ 6.8 x 10~ 3 . For such a value of v^l/yl) 
we take a minimum value of £, £ ~ 3.7 x 10~ 6 , under 
the perturbative unitary bound ffj/An < 1. Therefore, 
for a typical value = 500 GeV (£^,=0.033), together 
with k = kq = 0.45, we estimate £1 ~ 1.6 x 10~ 6 and 
£2^7.3xl0 -2 , which give the masses of the two charged 
scalars as Mi ~ 70 TeV and M2 — 335 GeV, respec- 
tively. As seen in Fig. 2, the choice n = kq = 0.45 gives 
\K/£\ max = 1-11 x 10 _1 , so that we predict 



Ren — Re 



Ru 



2.1 X 10 



-10 



(4.19) 



The values of Eq.(4.19) are too small to observe even 
in the near future colliders. Moreover, the case does not 
give any interesting neutrino phenomenology, because the 
neutrino mixing in this case comes out as 



U 



1 

"V2 






(4.20) 



\/2 V2 



with m Vl <C |m„ 2 | ^ \m V3 \. 

(B) feu ^ fre ^ fur'- The case is most interest- 
ing to us because it gives rise to bi-maximal mixing 
pattern as shown in Eq.(2.27) and in this case, the 
most serious bound on the /y coupling is given in 
Eq.(3.8). The relation among the fij coupling is given 
in Eq.(2.28) (which is necessary to explain reasonable 
value of Aml olar / Am1 tm ), so that it leads to the bounds 
as 



which predict 



V^fl/^l < 1.17 x 10- 2 , 

V^l/erl < 4.08 X 10- 5 , 

v ^|/ A1T |<2.84xlO- 7 , 



R ell < 5.2 x 10" 



(4.21) 
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R eT < 4.1 x 10 



-if) 



Rut < 1.8 X 1(T 14 , 



(4.22) 



where we have again utilized the constraints given in 
Eq.(4.21) giving £ ~ 1.1 x 10- 5 , \K/£\ max = 9.73 x 10~ 2 
of k = k = 0.44 and £ ^ 4.84 x 10" 6 (Mi ~ 41 TeV), 
ft ^ 7.5 x 10~ 2 (M 2 ~ 331 GeV). Unfortunately, simi- 
lar to the case (A), the values of Rij given in Eq.(4.22) 
are far away from the probing region of the future linear 
colliders. 

(C) / 2 T » / 2 ^ ^> / 2 e : In order to give an explanation 
of the excess value of /iAMM, we take VCl/url = 0.876 
from Eq.(3.20) and y/£\f Te \ ^ 5.3 x 1{T* from the 
Eq.(3.17) and y/£\f ei i\ = 4 -5 x 10_3 as obtained from 
Eq.(3.10). Utilizing the perturbative unitarity bound 
/ 2 T < 47T we obtain the value of £ as £ = 6.1 x 10~ 2 . For 
a typical choice of model parameters, M§ = 500 GeV, 
we obtain \K/£\ max ~ 9.46 x 10 -3 at k = k = 0.11, 
and we also get ft ~ 6.71 x 10~ 3 (M 2 ~ 1.1 TeV), 
ft ^ 3.02 x 10" 1 (M 2 ~ 166 GeV) . Therefore, we expect 
the maximal value of 



R eil < 1.6 x 10 



-12 



R eT < 1.1 x 10" 



R^r < 4.2 X 



10 



-16 



(4.23) 



Although we have obtained a large value only for R eT , 
the value R eT ~ 1.1 x 10~ 8 is still two order away 
from the observable region of the future linear collid- 
ers. Besides, since the case (C) gives the neutrino mixing 
given in Eq.(2.26) with the hierarchy of neutrino mass as 
\m Vl | <C \m V2 1 ~ \m V3 1, we must give up the explanation 
of the solar neutrino data within the framework of the 
three sequential neutrinos (y e , v^, v T ) in the Zee model. 

We summarized the main results of this section given 
in Eqs.(4.19), (4.22) and (4.23) in Table III. 



excess value of anomalous muon magnetic moment is in 
conflict with the hierarchical pattern / 2 M ^> / 2 e » / 2 T 
which is required to obtain bi-maximal neutrino mixing 
in the Zee model with three active neutrinos in order to 
explain the solar and atmospheric neutrino experimental 
results. We investigate LFV Z decays in the context of 
three different choices of hierarchical relations of the Zee 
fij coupling which are relevant to the present neutrino 
phenomenology. We find that among all the three decay 
modes of Z, only Z — > er decay gives rise to largest con- 
tribution to the ratio R eT ~ 10~ 8 which is two order less 
than the accessible value to be reached by the future lin- 
ear colliders only for the hierarchy of the /y coupling / 2 T 
S> ^> / 2 e as addressed in the case (C) which cannot 
reconcile the excess value of /iAMM and the solar neu- 
trino experimental result. Other possible LFV Z decay 
modes for all the three cases are significantly small to be 
observed in the next linear colliders. 

Although the values of the predicted branching ra- 
tios B(Z^e^eJ) are too small, the models (B) and 
(C) still remain as a promising candidate in the Zee 
model scenarios. The case (B) can give the simultane- 
ous explanation of the solar and atmospheric neutrino 
data, i.e., the bi-maximal mixing (2.26) together with 



solar 



/Am 



2 

atm 



V^me/m^ = 6.7 x 10~ 3 , (2.29), al- 
though it cannot give an explanation of the observed ex- 
cess of pAMM. The excess may be understood by the 
contributions from SUSY partners. On the other hand, 
the case (C) can give the simultaneous explanation of 
the atmospheric neutrino data and the observed excess 
of /iAMM, although it cannot give an explanation of the 
solar neutrino data. The solar neutrino data may be un- 
derstood H by an extended scenario with a sterile neu- 
trino |l(J, and so on. In the case (C), the neutrino mass 
matrix M v is given by the form 



E\ e 2 
M v = M„ T | ei 1 
£2 1 



(5.1) 



which leads to 



V. SUMMARY 



mi 



-2ei£ 2 M MT , 



In summary, we calculate LFV Z decays in the context of 
the Zee model and see their observability in the future lin- 
ear colliders subject to the existing bounds obtained from 
the \x decay, radiative /i and r LFV decays and /iAMM an( j 
experimental result. We first constrain the parameter 
space by considering the bounds obtained from the /i 
decay, so that the total contribution to the four Fermi 
coupling G p should not exceed the measured experimen- 
tal limits. We further have considered the bounds on 
the Yukawa couplings from the radiative charged lepton where 
decays and /iAMM. Both these constraints are given in 
Table I. The upshot of our analysis is that the hierar- 
chy of coupling / 2 T ^> / 2 M 3> / 2 e needed to explain the 



m 2 ~ -(1 - £i£ 2 )M MT , 



m 3 ~ (1 + £i£ 2 )M AI 



U 




f 



1 Ut \m T J 2 f qj,t 



(5.2) 



(5.3) 



T 2 - • (5.4) 
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If we give up to obtain a large value of B(Z^er), the 
parameter k becomes free from the constraint k = kq- 
Then, since M MT is given by 



AV ~2 



V6 



tan/3 m 2 . 



(5.5) 



we can obtain an arbitrary small value of l-M^I by taking 
k -> 1 (M 2 -* M ), with keeping the value V^/ MT =0.876. 

In conclusion, searches for LFV Z decay are not use- 
ful to confirm the two interesting models, (B) with 

fa » /re » fir and (C) with fa » fa » fa in the 

Zee's scenarios even at the near future colliders, where 
the model (B) can give the bi-maximal neutrino mixing 
together with the successful relation Am 2 oZar /Am 2 tm ~ 
y/2m e /m^ (although it fails to explain the observed ex- 
cess of /(AMM) and the model (C) can give the simulta- 
neous explanation of the atmospheric neutrino data and 
the excess of /iAMM (although it cannot give any expla- 
nation of the solar neutrino data within the framework of 
the three-flavor active neutrinos {v e , v^, v T )). Rather, an 
observation of the charged scalar (Zee scalar) associated 



with lepton flavor violation decay will be important in 
the near future colliders. 
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APPENDIX A 



Behavior of and 

The function F(6) is given by 

F(6) = / dx f dyx\n[l - ^xy(l - y)] 
Jo Jo 



1 



(4 - 56)6 + 16 tan" 



4- 6 



4V4^V / ^(2-6)tan- 1 



4-6 



(Al) 



For 6 -C 1, we obtain 



(AS) 



The function G(6) is given by 



GK ->=-H 



i . /-\_ ^ 2 -2(i-x) , i + V - ^ 



f/.c • — — — In- _ 

V6 V^ 2 - 4(1 - .*) 1 - yJ* iX 'lS~ xS - 



2(1 - .t) - ^a; 2 i 
+2 / dor ■ _ , ; s tan" 1 



V6 a/ 4 - 4a; - 6^ 2 



4(1 - a:) - 6a; 2 



(A3) 



where xo = (2/6)(vT+~6 ~ !)■ The approximate ex- 
pression of G(6) obtained as 

G(6) = -2.08276 + 3.06876 2 -2.97876 3 + 1.15536 4 --- • 

(A4) 



We can rewrite the function if(6> 6) explicitly as fol- 
lows: 



H(^,^) = j\x j\yx^\n 



1 - 



(1 - 6/6)y 
1 - 6a;y(i - y) 



+ ln 



1 - 



(1-6/6)3/ 
1 - 6a;y(l - y) 
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1 ,1 

dx / dy-xh\ 



ii 



(a-6) 2 2/(i-y) 



Ki - £16*2/(1 - 2/)] [6 - 66*2/(1 - y)] 



i ,1 

eta / dy-xh\ 
Jo 



l + B 



2/(1 - 2/) 



[l-&sy(l-y)][l- 6*2/(1 -2/)] 



(A5) 



where 



(6 - 6 
66 




(A6) 



For 6/6 < 10 2 , since £? ~ 6/6 > 10 2 , so that we can 
approximate i?(6;6) a s 

ff(6,6)-/ / dyx ln[l+Sy(l-y)]-F(^)-F(6) 
Jo Jo 

2 ^ V B Jl + I/B + l) 18 Va w 



.-l + iln|-i(£ 1 + £ 2 ). (AT) 

Note that the expression (A. 7) is valid only for 6/6 < 
10- 2 . 



[1] Y. Fukuda et al., Phys. Lett. B 335, 237 (1994); 
SuperKamiokande Collaboration, Y. Fukuda et al., 
Phys. Rev. Lett. 81, 1562 (1998); H. Sobel, Talk 
presented at Neutrino 2000, Su dbury, Canada, 2000, 



http:/ /nu2000.sno.laurentian.ca/ 



[2] Y. Suzuki, Ta lk presented at Neutrino 2000, S udbury, 



Canada, 2000, tittp://nu2000.sno.laurentian.ca/ 



[3] M. Cell-Mann, P. Ramond and R. Slansky, in : Su- 
pergravity, Proceedings of the Workshop, Stony Brook, 
New York, 1979, eds. P. Van Nieuwenhuizen and D. Z. 
Freedman (North Holland, Amsterdam, 1979) p. 315; T. 



[4] 



[5] 
[6] 

[7] 

[8] 
[9] 

[10] 

[11] 
[121 



[13] 



Yanagida, in: Proceedings of Workshop of Unified The- 
ory and Baryon Number in The Universe, eds. O. Sawada 
and A. Sugamoto (KEK 1979); R. N. Mohapatra and G. 
Senjanovic, Phys. Rev. Lett. 44, 912 (1980); Phys. Rev. 
D23, 165 (1981). 

C. D. Froggatt, M. Gibson and H. B. Nielsen, Phys. Lett. 

B 446, 256 (1999); R. Barbieri, L. J. Hall and A. Stru- 

mia, Phys. Lett. B 445, 407 (1999); R. Barbieri et al., J. 

High Energy Phys. 12, 017 (1998). 

A. Zee, Phys. Lett. 93B, 389 (1980). 

A. Yu. Smirnov and M. Tanimoto, Phys. Rev. D55, 1665 

(1997). 

C. Jarlskog, M. Matsuda, S. Skadhauge, M. Tanimoto, 
Phys. Lett. B449, 240 (1999). 

P. Frampton and S. L. Glashow, }iep-ph/9906375| (1999). 
S. Kanemura, T. Kasai, G.-L. Lin , Y. Okada, J.-J. Tseng 
and C.-P. Yuan, |hep-ph/0011357j (2000) . 
N. Gaur, A. Ghosal, E. Ma, P. Roy, et. al, Phys. Rev. 
D58, 071301 (1998). 

Y. Koide and A. Ghosal, Phys. Rev. D63, 037301 (2001) 
Muon (g-2) Collaboration: H. N. Brown, et al., friep- 
ex/0102017t Phys. Rev. Lett. 86, 2227 (2001). 
M. Anwar Mughal, M. Sadiq a nd K. Ahmed, Ph ys. Lett 
B417,87 (1998); E. J. Chun, hep-ph/0101170|, J. Bor- 
des, C Hong-Mo, T. S . Tsun, |hep-ph/0012lT9f ; E. O. II- 
tan, |iep-ph/0101017j ; S. Bergmann, H. V. Klapdor 



[14] 

[15] 
[16] 
[17] 

[18] 
[19] 



Kleingrothau s and H. Pas, |hcp-ph/0004048| ; J. I. Illana, 
T. Ri emann, |hep-ph/0010193 ; K. Cheung and O. C. W. 
Kong, ^ep-ph/0101347[ 

D. A. Dicus, H. J. He and J. N. Ng, |hcp-ph/0103126; ; E. 



Mituda and K. Sasaki, iep-ph/0103202 



A. Ghosal, Phys. Rev. D62, 092001 (2000). 

S. T. Petcov, Phys. Lett. B115, 401 (1982). 

D. E. Groom et al., (Particle Data Gr oup), Eur. Phys. 

Jour. C15, 1 (2000), |http://pdg.lbl.gov| . 

A. Czarnecki and W. J. Marciano, Nucl. Phys. (Proc. 

Suppl.) B76, 245 (1999); M. Davier and A. Hocker, Phys. 

Lett B435, 427 (1998). 

F. J. Yndurain, hep- ph/010231^ . 



10 



Table I 













fi — > eViVj 


1.37 x 1(T 4 


4 x Kn a 


4 x 1CT 3 


4 x 10~ ;l 


ei -> ej-y 




4.67 x 10 _& 


5.24 x 1(T 2 


3.39 x lCT^ 



TABLE I. Experimental upper bounds on £fikfkj {i,j,k = e,^i,r) from fi — > e^^j and — > e.,7 decays as discussed in 
Section II. A and B. 



Table II 







V?l/er| 






Case (A) : |/e/i.| — |/er| — 


0.68 x lO - * 


0.68 x 10^ 


0.68 x 10~ 2 


3.7 x 10~ b 


Case (B) : |/ eM | > |/ eT | > |/ MT 


1.17 x ltr 2 


4.08 x 10~ 5 


2.84 x 10~ Y 


1.1 x 10~ b 


Case (C) : |/ MT |» |/ e „| » f„\ 


4.57 x 10 - ' 1 


5.33 xl0~ ; ' 


0.876 


6.1 x W~* 



TABLE II. Possible maximal values of fij coupling under the constraints given in Table I and the result of /1AMM experi- 
ment. Also, the minimum value of £ under the constraints 7^ /471 < 1 is listed for each case. 



Table III 







max 


M 2 


Mi 




ReT 


RfJ,T 


Experimental upper bound 










5.05 x 10~ ; ' 


2.91 x 10~ 4 


3.56 x 10~ 4 


Case (A) 


0.45 


1.1 x 10 _I 


335 GeV 


70 TeV 


2.1 x 10~ iu 


2.1 x 10~ iu 


2.1 x 10~ lu 


Case (B) 


0.44 


9.7 x 1Q- Z 


331 GeV 


41 TeV 


5.2 x lO^ 4 


4.1 x 10~ iu 


1.8 x 10~ 14 


Case (C) 


0.11 


9.5X10" 3 


166 GeV 


1.1 TeV 


1.6 x 10~ ix 


1.1 x 10" 8 


4.2 x 10~ lb 



TABLE III. Prediction of the maximal values of Rij for the three different cases mentioned in Table II with the choices of 
model parameters as discussed in the text. The parameter values Mi and M2 are also listed. These value are fixed from the 
value of ko at which \K/£\ takes the maximum \K/^\ max . 
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FIG. 1. Z — *■ fi~r + decay in the Zee model due to charged scalar exchange. Z — > fi~e + and 
Z — > e~r + decays also arise due to the same diagrams with appropriate replacement of the internal 
fermion. 



10' 



10" 



M A = 1 TeV 




M d = 1 TeV 



M.= 200 GeV 



0.5 



K 



FIG. 2. The factor \K/£\ is plotted against k for typical 
cases (A) (solid lines), (B) (dotted lines) and (C) (dashed 
lines). 
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